Gravitational lensing of transient neutrino sources by black holes 
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Abstract 

r*j \ In this work we study gravitational lensing of neutrinos by Schwarzschild black holes. In 

Oh' particular, we analyze the case of a neutrino transient source associated with a gamma-ray 

burst lensed by a supermassive black hole located at the center of an interposed galaxy. We 
$H ' show that the primary and secondary images have an angular separation beyond the resolution 

C/3 | of forthcoming km-scale detectors, but the signals from each image have time delays between 

them that in most cases are longer than the typical duration of the intrinsic events. In this 
way, the signal from different images can be detected as separate events coming from the very 
£NJ ■ same location in the sky. This would render an event that otherwise might have had a low 

signal-to-noise ratio a clear detection, since the probability of a repetition of a signal from the 
same direction is negligible. The relativistic images are so faint and proximate that are beyond 
| the sensitivity and resolution of the next-generation instruments. 

^ ; PACS numbers: 95.30.Sf, 04.70.Bw, 98.62.Sb, 98.70.Rz. 

■ Keywords: gravitational lensing, black holes, neutrino astrophysics 

OO ' 
O ■ 

1 Introduction 

Gravitational lensing of photons by black holes has received great attention in the last few years, 
mainly due to the increasing evidence of the presence of supermassive black holes at the center of 
galaxies. Theoretical studies of black hole lenses, both numerically and analytically, were made 
with Schwarzschild [1-3], Reissner-Nordstrom [4], general spherically symmetric [5,6] and rotat- 
ing [7, 8] geometries, and also for black holes coming from alternative theories [9] or braneworld 
cosmologies [10]. Even naked singularities were considered as lenses [11]. Photons (or null mass 
particles) passing close enough to the photon sphere of the lens will have large deflection angles, and 
they can even make one or more turns around the deflector before reaching the observer. By this 
mechanism, two infinite sets of strong deflection images, one at each side of the lens, are produced. 
The presence of images with large deflection angles is not a new fact, since they were obtained 
already in 1959 for the Scharzschild spacetime [12]. The analytical study of these images is more 
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simple if one adopts the strong deflection limit, which consists in a logarithmic approximation of the 
deflection angle, first obtained for the Schwarzschild metric [12], revisited by several authors [2,13], 
extended to Reissner-Nordstrom geometry [4] , to general spherically symmetric spacetimes [5] and 
to Kerr metric [7]. For some lensing configurations two weak deflection images are also obtained, 
which are analyzed by making a first order Taylor expansion of the deflection angle (weak deflec- 
tion limit), as it is usually done for more standard astrophysical objects, such as stars and galaxies 
(see, e.g., [14]). Intermediate cases can be treated analytically by perturbative [15] or variational 
methods [16]. A special configuration, where no weak deflection images are present, is when the 
source is in front of the lens instead of behind it, which is called retrolensing [17]. Recently, the 
strong deflection limit was extended to include sources very close to the black hole [18]. 

Lensing of neutrinos have been previously studied by other authors. In Ref. [19], gravitational 
lensing of neutrinos by stars and galaxies was analyzed, and in Ref. [20], the lensing effects of 
supernova neutrinos by the Galactic center black hole was considered, in the weak deflection limit. 
However, perhaps the most interesting cosmological sources of neutrinos from the point of view of 
lensing are transients associated with gamma-ray bursts (GRBs). It is expected that proton-photon 
interactions during the GRB will result into copious photopion production and hence neutrinos 
would be generated from the decay of charged pions and muons (e.g. [21-24]). Since GRBs occur 
frequently, say once per day, and can be detected at gamma-rays by SWIFT satellite and then the 
follow up of the afterglows usually allows the identification of the host galaxy and the correspond- 
ing redshift (see, e.g., [25]), they are outstanding candidates for lensing produced by massive black 
holes in the center of interposed galaxies. 

We notice, however, that in the collapsar scenario for long GRBs [26, 27] the jet not always 
is expected to be able to make its way through the star, so no observable gamma-ray emission 
would result in such cases [22]. Nonetheless, the neutrino emission might be important. If the 
event is lensed, the neutrino signal should repeat and hence be identified, despite the absence of 
electromagnetic counterparts. 

In this letter we investigate gravitational lensing of neutrinos by Schwarzschild black holes. We 
pay special attention to neutrino transients lensed by supermassive black holes located at the center 
of galaxies. In Section [2] we present the expressions that give the positions and magnifications of 
the weak and strong deflection images, and in Section [3] we calculate the time delays between 
the arrival signals. Then, in Section 01 we calculate the specific time delays produced by some 
interposed supermassive black holes for neutrino transient at a distance of ~ 10 28 cm. Finally, in 
Section O a brief summary and the conclusions are presented. 

2 Positions and magnifications of the images 

Neutrinos have zero or negligible mass, so we assume that they follow null geodesies as photons 
do. We consider a point source of neutrinos, with angular diameter distance D os to the observer, 
behind a Schwarzschild black hole lens, placed at an angular diameter distance D Q \. The angular 
diameter distance between the lens and the source is dubbed D\ s . The optical axis is defined by 
the line that joins the observer with the deflector. The distances are very large compared to the 
Schwarzschild radius of the black hole and the angles are measured from the observer. We restrict 
our analysis to high alignment, which is more interesting from an astrophysical point of view, since 
the images are more prominent. Then the angular position of the source (3, taken positive here, is 
small. For this configuration, we have two weak deflection images and two infinite sets of strong 
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deflection (also called relativistic [1]) images. Neutrinos with closest approach distance ro much 
larger than the photon sphere radius r ps = 3MG/c 2 , which corresponds to the unstable circular 
orbit around the black hol^], will have a small deflection angle a, which can be approximated to 
first order in 1/tq by a = AGM/(c 2 ro) (weak deflection limit). Within this approximation, the lens 
equation has the form [14] 

where 9 is the angular position of the image and #e is the angular Einstein radius, given by 



with R s = 2MG/c 2 the Schwarzschild radius of the lens. The lens equation has two solutions: 



Pfi = -\B±y/fP + 4Bl), (3) 

that give the positions of the primary (upper sign) and the secondary (lower sign) images. The 
primary image lies inside the Einstein radius and the secondary image outside. When f3 = 0, 
instead of two images, an Einstein ring with radius #e is obtained. Another important aspect is the 
magnification of the images, defined as the ratio between the observed and intrinsic fluxes of the 
source. As a consequence of the Liouville theorem in curved spacetimes [31], gravitational lensing 
preserves surface brightness for neutrinos and photons, so the magnifications of the images are 
given by the ratio of the solid angles subtended by the images and the source, which result in [14]: 




where the plus sign corresponds to the primary image and the minus sign to the secondary one. If 
the position of the source (3 is close to zero, the magnifications of both images are large. If (5 = 
the approximation of point source breaks down and the magnifications become infinite. It is not 



difficult to see that fx p > 1 for all (3, and fi s > 1 only if (3/9 E < \] (3\/2 - 4)/2 » 0.35. When /3/0 E 
is large we have that /U p ~ 1 and fi s ~ 0. 

Besides the weak deflection images, two infinite sets of relativistic images are formed by neu- 
trinos that make one or more loops, in both directions of winding, around the black hole lens. For 
high alignment, the deflection angle corresponding to the relativistic images is close to an even 
number of it, a = ±(2mr + Aa n ) with < Aa n <C 1, the upper sign corresponding to one set 
of images and the lower one to the other set. The other angles involved are small, then the lens 
equation [ljl 

As 



tan (5 = tan 9 — — — (tan 9 + tan(a — 9)) , (5) 



D os 



takes the form [2,5] 



= 9 T ^Aa n . (6) 



For a complete study of the photon sphere in a spherically symmetric geometry see [28]; and for a general 
definition of the photon surface in an arbitrary spacetime see [29]. 

2 Eq. ((5]l is valid for asymptotically flat spacetimes, with the source and the observer in the flat region; for more 
general lens equations see [30]. 
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In the strong deflection limit, i.e. for trajectories passing close to the photon sphere of the black 
hole, the deflection angle can be approximated by a logarithmic function of the impact parameter 
b, defined as the perpendicular distance from the deflector to the asymptotic path at infinite. For 
the Schwarzschild geometry, it can be shown that [2, 5, 12] 



a 



-r-j In ( ^- - 1 ] + <■-, 



+ 0(b-b ps ), (7) 



with ci = 1, C2 = ln[216(7 — 4\/3)] — vr and 6 ps = V3r ps = 3\^3R s /2 the critical impact parameter. 
Neutrinos with impact parameter smaller than the critical value will spiral inside the photon sphere 
into the black hole, not reaching the observer, and those with b larger than b ps will make one or 
more outward turns outside the photon sphere, finally getting to the observer. As in the case 
of photons, using that b = sin D Q \9 9D Q \, inverting Eq. ((TJ) and Taylor expanding it around 
a = 2mr to obtain Aa n , then replacing the result in the lens equation (|6|) and finally inverting it, 
the positions of the relativistic images can be approximated (keeping only the lower order terms) 
by [2,5]: 

On = ±<£ + #7^/?, (8) 



where 
and 



(c2-2ri7r)/ci 

5 



bpS ^--^#^e n )(l + e n ), (9) 



it 



D i V AsA>ici 

is the n-th relativistic Einstein ring radius. For perfect alignment an infinite sequence of Einstein 
rings is obtained instead of point images. With the same considerations given above for the weak 
deflection images, the magnification of the n-th image has the same expression that was found 
previously for photons [2,5]: 



1 b 2 ps D, 
D^Dbd 



= -^TvTR (l + e n )e n , (10) 



for both sets of relativistic images. The first image (n = 1) is the strongest one and the others 
have magnifications that decrease exponentially with n. For a given source angle /3, the relativistic 
images are very faint compared with the weak deflection onea_j. 



3 Time delays 



Neutrinos that form distinct images take different paths, resulting in time delays between the 
images. Considering again that neutrinos follow null geodesies as photons do, the time delay 
between the primary and the secondary images is given by [14]: 



At 



p,s 



(1 + Zd) 

c 



2\9 P 6 S 



2 

^r+ln 



where Zd is the redshift of the deflector. The last equation can be written in the form 

2Rs 



At 



p,s 



1 + Zd) 



-pup 2 + m 



291 



+ ln 



P - ^(5 2 + 40. 



j3+Jf3' + A9i 



11) 



(12) 



3 For example, if /?/0e <C 1 we have that fAi/fip oc (Rs/D a i) 3 ^ 2 , which is usually a very small number. 
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When (3 = there is no time delay. Large time delays can be obtained if (3/0-e 3> 1, but in this case 
the magnification of the primary image is close to one and the secondary image is very faint. The 
optimal situation for a variable source is when (3/9e is small enough to have large magnifications 
of both images, but not too close to zero, so the time delay can be longer than the typical time 
scale of the transient source. 



In the case of relativistic images, the time delay between the images formed at the same side 
of the lens is given by [60: 



A^ m = -f(l + *0 



o i \ , o Ri \ _l V^ D os(w m - w n ) 
2-7r(n — m) + 2V2{w rn — w n ) ± — (3 



ClAs 



(13) 



where 



w k = e 



(c 2 -2&7r)/(2ci) 



and the upper/lower sign corresponds if both images are on the same/opposite side of the source. 
For the images at the opposite side of the lens we have [6]: 



At° jm = -^(l + z d ) 



/ \,o/o( \ I { y f 2D os (w m + w n ) 2D os 

2ir(n -m) + 2V2(w m - w n ) + — 

y ciAs As 







(14) 



where the image with winding number n is on the same side of the source and the other one on 
the opposite side. The first term in Eqs. (|13j) and (|14p is by large the most important one [6]. The 
time delays between the relativistic images are longer than the time delay between the primary 
and the secondary images. 



4 Lensing of neutrino transients 

The angular resolution of the primary and secondary images is beyond the capability of current 
and near future neutrino detectors, which is of the order of one tenth of a degree, but the temporal 
resolution of the images of individual transient events, which have typical durations in the range of 
~10 s to 100 s for long GRBs [25], is possible. As an example of this, we consider neutrino transients 
acting as possible sources situated at distances of the order of 10 28 cm, with supermassive black 
holes at the center of interposed galaxies as lenses. Some results of our calculations for specific cases 
of lenses in the local universe (z,i ~ 0) are shown in Tables [T] and [21 with the masses and distances 
taken from Ref. [32]. We see that the separation between the primary and secondary images of 
neutrino transients is of the order of a second of arc, so they cannot be resolved. For suitable values 
of the parameters involved, the weak deflection images can be both magnified several times, with 
time delays of 10 2 — 10 4 s, larger than the intrinsic time of variation of the sources. 

If one fixes (3 /#e to obtain from Eq. @ the desired values of magnification of the images, it 
is clear from Eq. (|12f) that the time delay increases linearly with the redshift of the lens. Then, 
with a typical source with redshift z s ~ 1, the values of time delays between the primary and the 
secondary images for lenses closer to the neutrino sources can be up to twice of those obtained in 
Tables [1] and [2j But far lenses require better alignment to have large magnifications because #e 

4 The expressions from Ref. [6] have been rewritten here using physical units, adding the cosmological factor 1 + Zd 
and expanding them to first order in the source position angle (measured from the observer instead of from the 
source) . 
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Galaxy 


Black hole mass (Mq) 


Distance (Mpc) 


9e (arcsec) 


Ats, P (s) 


Milky Way 


2.8 x 10 6 


0.0085 


1.6 


11 


NGC0224 


m 

3.0 x 10 7 


0.7 


0.6 


1.2 x 10 2 


NGC3115 


2.0 x 10 y 


8.4 


1.4 


7.9 x 10 3 


NGC3377 


1.8 x 10 s 


9.9 


0.4 


7.1 x 10 2 


NGC4486B 


5.7 x 10 8 


15.3 


0.5 


75 — 

2.2 x 10 3 


NGC4486 


3.3 x 10 y 


15.3 


1.3 


1.3 x 10 4 


NGC4261 


4.5 x 10 s 


27.4 


0.4 


1.8 x 10 3 


NGC7052 


3.3 x 10 s 


58.7 


0.2 


1.3 x 10 3 



Table 1: Time delays between the weak deflection images of neutrino burst sources at a distance 
of 10 28 cm. The lenses are supermassive black holes at the center of the galaxies indicated. The 
Schwarzschild geometry was adopted to model the black holes. The source angular position is 
= 0.1 #E; with #e the angular Einstein radius. In this case, the angular positions of the primary 
and the secondary images are 9 V = 1.05 #e and 9 S = — 0.95 #e, while their respective magnifications 
are \x p = 5.5 and ^ s = 4.5. 

decreases with the distance to the lens. 

Throughout this letter we have assumed that the neutrinos move like photons in null geodesies, 
so the time delays do not depend on the energy of the neutrinos. The time lag between neutri- 
nos with an energy E v and a small mass m u and photons travelling the same distance d can be 
approximated by At {1 /2){d / c){m y c 2 / E u ) 2 , which using E v > 1 TeV (for neutrinos associated 
with GRBs), m u c 2 < 1 eV, d ~ 10 28 cm, gives At < 10~ 6 s. This time lag corresponds to a total 
travelling time of about 3 x 10 17 s. Then, if neutrinos have mass, the time delays given in Tables 
[T] and should be modified in the same proportion, i.e. in less than one part in 10 23 . There is 
also observational evidence related to the supernova SN1987A which shows that the time delay 
due to the presence of the galaxy for photons and neutrinos with different energies is the same 
within a 0.5% or better accuracy [34]. All these justify our assumption, and the results obtained 
are excellent approximations if the neutrinos have sub-eV mass. 

Concerning the relativistic images, if we choose 0/9e = 0.1 as it was done in Table [TJ we have 
from Eq. (jlOp that the magnification of the brightest image is [i\ = 1.1 x 10 -17 for the Galactic 
black hole and = 1.8 x 10 -22 for NGC4486. Similar values are obtained for the other black 
holes considered in Tables [T] and [5J The other image magnifications decrease exponentially with 
n. To obtain magnifications about one or larger, a closer alignment is necessary, i.e. ~ b ps /d \ 
instead of ~ #e- Then, while the primary and secondary images are amplified several times for 
0/&K = 0.1, the strong deflection ones are highly demagnified. Using Eq. (jHJ), it can be seen that 
the angular separation between the strong deflection images is of the order of micro arc seconds or 
less. The time delays between the relativistic images, given by Eqs. (|13p and (|14f) . can be large, but 
they are too faint to be detected. So, in what follows, we restrict ourselves to the weak deflection 
images. 

The probability of supermassive black holes located at the center of galaxies in the line of sight 
to GRBs is not negligible because of the high-redshift of most GRBs. Moreover, optical spec- 
troscopic observations can detect absorbing lines of the interposed galaxy in the afterglow, hence 
allowing a direct determination of the different distances involved in the scenario. In the case of 
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Galaxy 


Black hole mass (Mq) 


Distance (Mpc) 


9e (arcsec) 


At s ,p (s) 


Milky Way 


2.8 x 10 6 


0.0085 


1.6 


55 


NGC0224 


m 

3.0 x 10' 


0.7 


0.6 


— „ „ m — 

6.0 x 10 2 


NGC3115 


2.0 x 10 y 


8.4 


1.4 


4.0 x 10 4 


NGC3377 


1.8 x 10 s 


9.9 


0.4 


3.6 x 10 3 


NGC4486B 


5.7 x 10 8 


15.3 


0.5 


1.1 x 10 4 


NGC4486 


3.3 x 10 y 


15.3 


1.3 


6.6 x 10 4 


NGC4261 


4.5 x 10 s 


27.4 


0.4 


8.9 x 10 3 


NGC7052 


3.3 x 10 s 


58.7 


0.2 


6.6 x 10 3 



Table 2: Time delays between the weak deflection images of neutrino burst sources at a distance 
of 10 28 cm. The lenses are supermassive black holes at the center of the galaxies indicated. The 
Schwarzschild geometry was adopted to model the black holes. The source angular position is 
(3 = 0.5 with #e the angular Einstein radius. In this case, the angular positions of the primary 
and the secondary images are 9 V = 1.28 and 8 S = — 0.78#e> while their respective magnifications 
are // p = 1.6 and fj, s = 0.6. 

choked GRBs, where only neutrinos are produced, the time delays and the relative magnifications 
of the signals could used for the unequivocal identifications of dark neutrino transients. A neutrino 
transient associated with a GRBs might have an fluence of several times 10 -4 erg cm -2 [23]. With 
a mild amplification as obtained for the parameters adopted in Tables Q] and El this might imply 
the detection of a few neutrinos by a km 3 -detector. Even if the signal-to-noise ratio is not at a high 
confidence level, the repetition of the signal on a time scale from minutes to hours from the same 
location in the sky would render the identification of the neutrino transient source unequivocal. 
If the detection of the GRB afterglow allows a clear determination of the redshifts involved, then 
Eqs. Q and (|12|) can be used to obtain an independent estimated of the central black hole mass 
in the interposed galaxy. 

The analysis of current databases indicates that the space-time clustering of GRBs is only 
marginal, at the level of 5% or less [33], but as we have mentioned in the Introduction, choked 
collapsars can result in transient neutrino sources without electromagnetic counterparts, so the 
total number of neutrino transients that is affected by lensing effects could be significantly larger 
from what is inferred from GRB population studies. The detection of a single event could be of 
paramount importance for our understanding of physical processes governing the GRBs. 

5 Final remarks 

In this letter we have shown that the primary and secondary images of neutrino transient sources 
lensed by supermassive black holes cannot be angularly resolved but they could be temporally 
resolved by next generation instruments. The relativistic images, instead, are too faint and packed 
to be detected. Thus, we have found that neutrino transients produced by long GRBs can act as 
sources for gravitational lensing when supermassive black holes are present in foreground galaxies. 
This sources would have a unique signature, that will allow an easy detection above the background 
despite a possible low signal-to-noise ratio: repetition. The neutrino fluence can be magnified, but 
more importantly, the arriving signal will repeat, leading to an unequivocal identification. We 
conclude that neutrino gravitational lensing can help to establish GRBs as sources of relativistic 
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protons and neutrinos, as proposed by several authors [22,23]. 
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